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SUMMARY 

With the advent of lunar landings planned for the near future, it 
is essential that considerable though be given to the behavior of jet ex- 
haust in a vacuum under low gravity conditions, especially in relation 
to the most probable lunar surface conditions. Of primary concern is the 
motion of lunar surface particles resulting from the impingement of ex- 
haust gases from a lunar landing or takeoff space vehicle and the result- 
ing hazards to the vehicle itself and to both local and distant manned 
lunar bases or facilities. The combined parameters of low gravity, high 
vacuum, and high energy exhaust gases result in the displacement of large 
masses over considerable distances, so that the problem is not merely a 
localized one. 

This preliminary study was directed specifically toward determination 
of the trajectories or flight paths of lunar surface particles accelerated 
by the jet blast from a lunar landing or return vehicle. A range of par- 
ticle velocities was assumed based on an anticipated effective exhaust gas 
velocity of 13,250 ft/sec for a LOX-LH 2 engine. Thus, velocities both 
above and below those required for escape or orbit from the lunar surface 
were considered. Data are presented in graph and chart form, along with 
a brief discussion. Sample calculations and basic data for the ballistics 
study are presented in Appendix I. Appendix II describes in detail the 
follow-on work to be accomplished and the problems to be encountered in 
an analysis aimed at predicting the mass and distribution of particles 
most likely to be associated with the velocity range assumed in Appendix 
I. Determination of the momentum energy relationships between 

exhaust gases and lunar surface particles of assumed sizes and density 
will give us a better feel for the degree of hazard involved. 

As discussed later in Section II and Appendix II, large dust clouds, 
and the ballistics of relatively large lunar surface particles resulting 
from a lunar landing or return launching may present serious problems in 
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relation to the landing operation itself and most definitely can produce 
hazardous conditions for both distant and surrounding lunar base facilities, 
equipment, and personnel. 

Of course, more studies will be required in these areas when more 
definitive data on the exact composition and structure of the lunar surface 
is obtained. Such studies will be needed to provide data upon which to 
base lunar facility design concepts and operational criteria. 


SECTION I. DISCUSSION OF RESULTS 

It was considered reasonable to assume that lunar surface particles 
would attain velocities ranging from 0 to 8000 ft/sec when subjected to 
the exhaust of a LOX-H 2 engine with an effective exit velocity of 13,250 
ft/sec, considering that some velocity energy will be lost in transmission 
through standoff shock waves and through deflection from the lunar surface. 
Since the escape velocity from the lunar surface is 7800 ft/sec and the 
circular orbital velocity is 5540 ft/sec, one can readily envision that 
escape, orbiting of particles, or extremely eccentric partial orbits may 
occur as a result of exhaust jet impingement. Appendix I contains sample 
calculations and a summary of data (Tables I, II, and III) for a number 
of assumed initial particle velocities and initial angles with respect to 
the local vertical. The varying angles were considered because of assumed 
erosion of the lunar surface with time during jet impingement. Test data 
in Reference 1 includes that this pattern will be in the form of a conic 
section, with a final semi-vertex angle of approximately 50 degrees. 

Figure 1 shows the hyperbolic flight path for a particle accelerated 
to a velocity greater than escape. For assumed velocities of 7410 ft/sec 
and 6350 ft/sec, which are less than escape velocity, Figures 1 and 2 
apply. The velocities exceed circular orbital velocity at the lunar sur- 
face so that an eccentric orbit is possible. However, here we note the 
stringent requirements for an orbit initiated from the lunar surface. 
Assuming the particle has an initial velocity exceeding the circular 
orbital velocity, then the initial launch point would be the low point 
or perigee of a hypothetical orbit. Also since it is launched from the 
surface instead of at some altitude, the takeoff angle with respect to the 
local vertical would necessarily be exactly 90 degrees in order that an 
orbit be possible. The possibility of an orbit is further restricted in 
that the initial particle impulse would necessarily have to take place 
from a very high point so that the particle would not impact on a nearby 
hillside or other obstruction prior to traveling over the horizon and 
attaining or orbit. Thus, with a firing angle of 90 degrees and a ve- 
locity greater than circular velocity, one complete orbit is possible 
but highly improbable. At any angle other than 90 degrees, eccentric 
partial orbits resulting in long surface ranges would be experienced. 
Figures 1 and 2 describe the flight paths for initial velocities of 7410 
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and 6350 ft /sec and number of assumed angles with respect to the local 
vertical. Figure 3 depicts the actual lunar surface range covered by 
these same particles. Of course, for the departure angle of zero degrees 
from the horizontal, the theoretical distance given is equivalent to one 
orbit or, more precisely, the circumference of the moon. 

For assumed velocities less than circular orbital velocity, figure 
4 applies. In this case, conic section central force field ballistic 
equations taking into account the moon s curvature were utilized in de- 
termining the resulting particle ranges (Appendix I). Figure 5 describes 
particle range based on velocities from 88 to 550 ft/sec and as defined 
by standard ballistic equations neglecting curvature. 

Preliminary conclusions arrived from following this analysis will be 
discussed in the following section. 

SECTION II. RECOMMENDATIONS AND CONCLUSIONS 

A review of the data indicating the ranges attained by objects as 
a function of initial velocity leads one to conclude that depending upon 
surface composition, a lunar landing or return operation may result in 
the bombardment of both local and d istant lunar bases with relatively 
large size particles. This is especially evident under the existing 
lunar conditions of vacuum, low gravity, and high energy exhaust gases. 

Such bombardment could prove hazardous to facilities, associated 
equipment, and personnel, especially where lightweight inflatable thin 
film strucutes or shelters are proposed. These problems must be con- 
sidered in the siting and design of future lunar facilities so that ade- 
quate protection will be provided. 

Jet impingement may also result in the formation of large dust clouds 
covering a considerable area of the lunar surface, thus proving hazardous 
to the landing operation. Coverage of large surface areas is expected 
becuase of jet expansion effects resulting from the gross under -expansion 
of jet nozzles in a vacuum. As shown from data of a preliminary test 
data in a vacuum (Appendix II), cratering will occur upon jet impingement 
on a "soft” lunar surface. If the scale model tests conducted by Lewis 
are any indication, the ’’cratering' 1 from a full scale LOX-LH 2 nozzle will 
be quite sizeable and will present a formidable design problem in develop- 
ing a safe lunar landing system. Of course, all of this is based on 
assumed lunar surface conditions which may or may not exist. Initial 
unmanned landings will probably clarify this aspect of the problem. 

Because of the critical nature of these problems, it is concluded 
that additional study work should be directed toward determination of 
the particle masses that will most probably be accelerated to the range 
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of velocities just studied. This, of course, would necessarily be based 
on a range of assumed particle densities. The general approach to the 
future study work to be performed is discussed in Appendix II. It will 
first be necessary to determine the overpressures at the deflecting 
surface as a function of nozzle height above the surface, propellant 
type, and nozzle characteristics, such as operating chamber pressure, 
and available gas stream momentum energy must be related in parametric 
equations predicting probable particle motion. For both phases of the 
problem, and especially Part 2 (as discussed more fully in Appendix II) 
a considerable amount of additional test work is required. This is 
necessary because the variable particle velocity factor cannot be sepa- 
rated in an equation from the other time-dependent variables without 
making a number of assumptions based on the existing supply of inade- 
quate test data. 
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APPENDIX I - SAMPLE CALCULATIONS 
BASIC LUNAR DATA 


R = Radius of Moon = 1082 miles 

W = Angular rotation = 1 rev/27.33 earth days 

Surface Velocity at Equator = 10.4 miles/hr 

g = 1/6 g ^ = 5.37 ft/sec 2 

6 moon to earth 

M = Mass 
m 

Mass 


moon " 81.375 MasS earth 
= 6.595 x 10^ tons 


earth 

G = Universal gravitational force constant = 

PROPULSION DATA 


3.44 x 10~ 8 ft 4 
lb-sec^ 


Assume that the orbital launch or landing vehicle is propelled by two 
LOX-LH 2 engines, rated at 15,000 lb thrust each. 

Assume I or , of these engines in vacuum = 412 sec 
sp 

/ Thrust 15,000 ^ „ , 

Propellant flow/engine = — = — =* 36.5 Ib/sec 

I sp 

For both engines, total gas flow = 73 lb force/sec 

Effective Exhaust Velocity = I (g .,) 

J sp °eartlr 

= 412 sec (32.2 ft/sec 2 ) = 13, 250 ft/sec 

Based on this exhaust velocity, potential escape, orbit, and bal- 
listic flight paths of lunar surface particles must be considered. A 
sample calculation for various assumed velocities will now be presented. 


Case 1 - V. . . . . > escape velocity 
initial r J 


V from moon 

escape 


GM 


G = universal gravitational constant 
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M = Mass of moon 

r = distance from center of moon to point of injection = 

V = 7800 ft/sec at the lunar surface 

escape 

Since the assumed particle velocity exceeds the escape velocity at 
the lunar surface, it will escape from the moon and fall toward the earth 
for any angle of B 0 sufficient to clear local obstacles (Fig. 1). 

Case 2 - Assume velocity = 6350 ft/sec 

The initial velocity of this particle, 6350 ft/sec, is less than 
escape velocity but more than the velocity required for a circular orbit 
at the lunar surface. 


V 


circular orbit 


g R 


2 

moon 


orbit 


5540 ft/sec 


Since the initial particle velocity does exceed the orbital velocity, 
an orbit is potential with the low point or perigee being the point of 
takeoff. A potential orbit can also be verified by reference to the equa- 
tion defining X on Figure 6. 

V 2 r V 2 

„ O O O ? 

X = — — = where GM = g r 

GM g r & o 

° o 

X is a parameter equivalent to twice the ratio of kinetic to potential 
energy of the particle at takeoff and is also equal to 


Therefore, if X is = > 1, the local 
velocity is higher than circular 
orbital velocity, and an orbit is 
potential. 


actual 


\ 


reqd for 
circular orbit 


For the case of V = 6350 ft/sec, 

_ (6 . 35) 2 x 10 6 ft 2 sec 2 

sec 2 5.37 ft (1.082 x 10 3 ) (5.28 x 10 3 ) mile 


1.315 


7 


.*. , satellite orbit is possible, depending on angle of injection. 

If angle B 0 is 90 degrees, the eccentricity of a potential orbit is 
defined by 


= V 1 - X (2 - X) sin 2 


By, 


e = V 1 - 1.315 (2 - 1.315) (1) = 0.316 


Flight paths are defined by the following (Fig. 7): 


e < 1 - elliptical 
e = o - circular 


e = 1 - parabolic 

e > 1 - hyperbolic 

.*. elliptical equations apply. 

1 + e 

The apogee of the potential orbit is defined by r = (r Q ) 

max . z - A 

where r 0 in this case = radius of injection from lunar centroid = R moon 

1316 

apogee or r = ‘ QC (1082) = 2080 miles, or 1000 miles above the 

K & max. 0.685 ' 

lunar surface. 


Radius of perigee 


1 - e 

2 - X 


(r ) = 1 (R ) 
o moon 


= R 


moon 


Therefore, if B 0 is 90 degrees and the particle is initially launched 
from a high point on the lunar surface, it could conceivably make one 
complete orbit before impact. As could be expected, however, the possi- 
bility of an orbit is practically nil because of the stringent requirements. 

It can be shown that for angles of B Q less than 90 degrees calcu- 
lated perigee point is less than the radius of the moon so that no orbit 
is possible. To determine the actual range covered on the lunar surface 
(Fig. 7). 

r, which is shown as the distance to a point on the elliptical path 
is defined by the equation 


r = ^ 

1 + e cos a 
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p is the semilatus rectum of the ellipse 



a 


V ^ sin^ B 
o o 


g 


o 


a(l-e 2 ) 


a is the angle measured from the perigee to a point on the elliptical 
path (Fig. 7). 


The value of r as a function of e,p, and cos a then defines the 
elliptical path. In considering potential orbits, by letting r = R of 
the moon, angle a indicates the angle from the perigee to a point on the 
elliptical path a distance equal to the moon's radius from the foci. To 
have an orbit, angle a must be 0 degrees. Any angle of a greater than 
0 degrees indicates that the perigee radius of the theoretical elliptical 
orbit is less than the radius of the moon, which, of course, makes a 
closed orbit impossible. For example, consider the case for V q = 7410 

ft/sec. By firing horizontally, (B = 90 degrees), P = 1940 miles, and 
e = 0. 79, • 

P 

r = 

1 + e cos a 


where 


r = R 


moon" 


1 

cos a = — 
e 



.*. a = 0° (or 360°) 


This indicates that the point of particle take-off from the lunar 
surface is the perigee of the orbit, and an orbit could occur if the 
particle did not hit any local obstructions prior to traveling over the 
local horizon (Fig. 1). 


It can be shown that for any angle of B q less than 90 degrees, a 

complete orbit is not possible. Consider V - 7410 ft/sec with B = 30 
degrees. 2 2 

b 2 V q sin B q 2 


p = - = 


g 


= XZ^_x_10li0^5I = 2> 56 x 10 6 ft 


e = eccentric 


ity = -J. 


1 - X (2 - X) sin 30° = 0. 952 


let r = R = 5. 7 x 10 ft 
moon 


r = f 


1 + e cos a 
a = cos 


1 [i <£ - 1)] for [ ] > 0 
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a = jt - cos 
a = 125. 45° 



> 0 


Thus, no orbit is possible. Since the point of departure from the 
lunar surface remains constant, this indicates that the hypothetical 
perigee of the ellipse and the major axis of the ellipse has shifted an 
angle a of 125.45 degrees with respect to the ellipse where B = 90 
degrees. (Fig. 1, axis shift from x-x to z-z). Since the ellipse is 
symmetrical, the impact point on the opposite side of the axis will also 
be defined by angle ql> and the angular range covered on the lunar surface 
by the partial orbit will be defined by 360 degrees - 2 OL. 


For this case, the surface range(s) is defined as follows: 


(s) = R (d 0) = R ^ 36 ° ' 2 ir 

moon moon 180 


(s) 


1082 miles 


36 0-2 (125.45) 
180 


it = 2060 miles. 


After determing the axis shift for each angle of departure (B) , the 
relationships established above were utilized in plotting the flight 
paths shown in Fig. 1 and 2 for velocities of 7410 and 6350 ft/sec 
(See Fig. 3 also). 

Case 3 - Assume V = 4040 ft/sec (Fig. 4) 

An initial velocity of 4040 ft/sec is less than escape or orbital 
velocity. Also, X is considerably less than 1. Therefore, ballistic 
range equations are applicable. As given by Fried & Richardson (Ref. 1) , 
the range angle 0 (the circular angle transversed by the object from 
launch to impact) is given by the following equation; 



2 


r , . .2 i 

. -1 
sm 

(1-crXsinB) i 

o 

-1 

- Tan 

1-X sin B 

o 

(1 - 2 X sin^ B + X^ sin^ B )4 
v o o _ 

sin B cos B 


where o = 1 + — 


*h is defined as the altitude at burnout. In our case, which is identical 
to impulsive burning or instantaneous application of all thrust at the 
lunar surface, h = 0. /. a = 1 
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The range from take-off to impact may also be defined as follows: 


Range = 2 R tan -1 

0 mnnn 


moon 


V L sin B cos B 
o o o 


L R 


moon moon 


- V 2 sin 2 B 


For B = 75 c 
o 


Range = 5.7 x 10° ft (2) tan 


-1 


(4.04)^ x 10 6 (0.9659) (0.2588) 

L 30.6 x 10 6 - (4. 04) 2 x 10 6 (0.934 J 


tan 


_1 = 0.266 = 14.9° 


ft 


. _ 11.4 x 10^ mile 14.9 it 

Range ^0 = Similes 


A similar approach was taken for other angles of B and for reduced 
velocities. When the range became very short (V < 1000 ft/sec) the use 
of basic ballistic equations neglecting curvature was found acceptable 
(Fig- 5). 


PERI OP OF REVOLUTION -ELLIPTICAL ORBITS 


The period of revolution (Figs. 1 and 2) of the hypothetical particle 
orbits where the initial velocity exceeded orbital velocity is defined by 
the following equation: 

where a = - — ^ 7 — (Fig. 7) 

and [i = gravitational factor = gR^ 

o 

g = gravitational pull - ft/sec^ 

R = radius from center of lunar mass to orbit 

Data for all of the above referenced figures is summarized in 
Tables I, II, and III. 



TABLE I. LUNAR BALLISTICS DATA 


II 





TABLE II. RANGE VERSUS ANGLE OF DEPARTURE DATA 
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60° 

560 

268 

102 

44 

24.8 

11.6 

11 

9.25 

7.64 

4.89 

2.75 

1.22 

0.306 

0.237 


• 

O 

m 

750 

332 

120 

51 

28.6 

13.4 

12.7 

10.7 
8.8 
5.65 
3.18 
1.41 
0.353 
0.274 

r — 1 

Ctf 

u 
a) a 
• u o 

W 3 N 

o 

O 

CO 

795 

313 

109 

44 

24.8 

11.6 

11 

9.25 

7.64 

4.89 

2.75 

1.22 

0.306 

0.237 

(Miles) \ 
of Depart 
Local Hori 

0 

m 

i — i 

560 

194 

64 

25.4 

14.3 

6.7 

6.35 

5.35 
4.4 
2.82 
1.59 
0.7 
0.174 
0.137 


Range 

Angle 

with 

o 

O 

o 

ii 

Init. Velocity 
(Ft/sec) 

v 


V = 4040 
o 

2870 

1800 

1200 

900 

616 

600 

550 

500 

400 

300 

200 

100 

88 

1 
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TABLE III. FLIGHT PATH DATA (V Q > CIRCULAR ORBITAL VELOCITY) 
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TABLE III. FLIGHT PATH DATA (V Q > CIRCULAR ORBITAL VELOCITY) (Cont ' d) 


V Q = 7410 ft/sec 


B_(Degrees) a (Degrees)l e_ |_ 2 _(Hiles)| r_(Miles)| e cos al cos 


(Axis Shift 

66 ) 


-0.847 


+0.847 


0.900 970 


0.845 


0.636 


-0.636 


-0.636 


+0.636 


30 

(Axis Shift 
125.45) 


0.952 485 


0.895 


0.673 


-0.673 


10,000 


-0.952 


-0.673 
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TABLE III- FLIGHT PATH DATA (V > CIRCULAR ORBITAL VELOCITY) (Cont'J) 

= 7410 ft/sec 


B (Degrees) 

q (Degrees) 

e 

p (Miles) 

r (Miles) 

e cos a 


270 


485 

485 

0 


315 


485 

290 

+0.673 


360 


485 

248 

+0.952 
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TABLE III- FLIGHT PATH DATA (V Q > CIRCULAR ORBITAL VELOCITY) (Cont'd) 

V = 6350 Ft/sec 
o 


B (Degrees) 

q (Degrees) 

e 

2 _(Miles) 

r (Miles) 

e cos a 

1 

90 

20 

0.316 

1420 

1095 

0.297 



45 


1420 

1160 

0.224 



90 


1420 

1420 

0 



135 


1420 

1830 

-0.224 



180 


1420 

• 

2080 

-0.316 



225 


1420 

1830 

-0.224 



270 


1420 

1420 

0 



315 

0.316 

1420 

1160 

+0.224 



360 


1420 

1080 

+0.316 


75 

(Axis Shift 

20 

0.400 

1330 

966 

0.376 


55) 

45 


1330 

1035 

0.283 



90 


1330 

1330 

0 



135 

i 

1330 

1855 

-0.283 



180 


1330 

| 

2220 

-0.4 



225 


1330 

1855 

-0.283 



270 


1330 

1330 

0 



315 


1330 

1035 

+0.283 



360 


1330 

950 

+0.4 


60 

20 

5.700 

1065 

695 

0.535 


(Axis Shift 

45 


1065 

759 

0.403 


91.2) 

90 


1065 

1065 

0 



135 


1065 

1780 

-0.403 
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TABLE III- FLIGHT PATH DATA (V Q > CIRCULAR ORBITAL VELOCITY) (Cont’d) 

= 6350 ft/sec (Cont 1 d) 


B (Degrees) 

a (Degrees) 

e 

p (Miles) 

r (Miles) 

e cos a 








180 


1065 

2480 

-0.57 


225 


1065 

1780 

-0.403 


270 


1065 

1065 

0 


315 


1065 

759 

+0.403 

* 

360 


1065 

679 

+0.57 

45 

20 

0.742 

713 

420 

0.697 

(Axis Shift 

45 


713 

467 

0.525 

117.25) 

90 


713 

713 

0 


135 


713 

1500 

-0.525 


180 


713 

2760 

-0.742 


225 


713 

1500 

-0.525 


270 

0.742 

713 

713 

0 


315 


713 

467 

+0.525 


360 


713 

409 

+0.742 

30 

(Axis Shift 

20 

0.880 

356 

195 

0.826 

139.55) 

45 


356 

220 

0.622 


90 


356 

356 

0 


135 


356 

941 

-0.622 


180 


356 

2960 

-0.88 


225 


356 

941 

-0.622 


270 


356 

356 

0 


315 


356 

220 

+0.622 

J 

360 

j 


356 

189 

+0.88 
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APPENDIX II 

PROBLEM AREAS FOR FUTURE STUDY 

This preliminary study has indicated the velocity ranges and 
ballistics to be encountered on blast accelerated lunar particles. The 
problem now remains to determine the masses of particles that will most 
probably be accelerated to these velocities. The problem of establishing 
the momentum energy transferred to these particles is fraught with many 
uncertanties and inseparable variables. With knowledge of exhaust gas 
densities, mass flow rates, stand-off shock wave theory, and gas velocity, 
one can determine the maximum over-pressure and resulting forces that could 
conceivable be initially applied to an object in the jet stream. Through 
knowledge of the initial temperature and flow rate of gases, the energy 
loss at impact, and the probable angles of deflection as a result of 
surface erosion, one can use a molecular approach to calculate the momentum 
energy within the gases per unit area of dispersion for a given distance 
or time and thus calculate the force that would be applied to a stationary 
object at that point. However, since the particle is in motion, the force, 
and of course the velocity and acceleration curves of the particle will 
be an integral function of time. In typical dimensional analysis, blast 
wave variables are determined as a function of time by four major para- 
meters: (1) Shock over-pressure, (2) Ambient pressure, (3) Time duration 

of the positive winds and (4) The speed of sound in the undisturbed air. 

Of course, items (2) and (4) become zero under lunar conditions. Further- 
more, it is soon realized that even for conditions such as that on earth, 
a completely pure mathematical solution is not possible because the vari- 
able particle velocity factor cannot be separated in an equation from 
the other time -dependent variables. For this reason, most blast effect 
equations here on earth are based on empirical relations determined by 
fitting curves to computed blast data (Ref. 2). 

Practically all existing theory developed on secondary-missile 
behavior has come about as a result of weapons effect studies. Very 
often, semi -empirical relations are used to define over-pressure dura- 
tion as a function of blast yield, over-pressure, ambient pressure, and 
the speed of sound. Then, with knowledge of the velocity of propagation 
of the pressure disturbance, the displaced object motion can be predicted. 
None of these analyses are adaptable to our problem, since there is not 
an instantaneous release of energy and there really is no wave or over- 
pressure propagation in a vacuum; there is only gas dispersion, the 
initial dispersion being a function of nozzle height above the lunar 
surface and the nozzle configuration or degree of under-expansion. 

A review of nearly all of the listed references relating to blast or 
jet exhaust effects indicates that most work to date is based specifically 
on one type of rocket engine and one specific situation, and is not trans- 
posable to another engine or operating condition. This indication is also 
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brought out clearly in Ref- 3. What is really required in the follow-on 
study to this report is a breakdown of the problem into the following 
areas : 


1. A determination of the over-pressure at the deflecting surface 
as a function of nozzle height above the surface, propellant type, and 
nozzle characteristics such as operating pressure, nozzle contour, and 
expansion ratio. 

2. Definition of the relationships between surface over-pressures 
and momentum energy transfer to secondary particles. Because of the diffi- 
culties involved in a pure mathematical approach to Part 2, it appears 
that it will be necessary to implement basic parameteric equations with 
some test data and/or empirical relationships. Perhaps the first effort 

at obtaining such data was made by Lewis Research Center (Ref. 4 ) where 
recent tests were made in a vacuum to determine the pressure distributions 
and erosion patterns on simulated hard and soft lunar surfaces. These 
model tests were performed with cold air jets at nozzle pressure ratios 
up to 288,000. Some of the more important results of these tests were: 

a. Surface pressure distributions were dependent on nozzle 

area ratio. For nozzles with area ratios ranging from 1 to 25 , the nozzle 
total pressure ratio remained constant at 288,000- However, the nozzle 
static pressure ratio P e xit/ p ambient > which is an index of jet spreading, 
varied from 152,000 for a sonic nozzle (A exit = A t ^ roat ) to 545 for a 
conical nozzle with an area ratio of 25 and a corresponding exit Mach 
number of 5. At the largest descent heights (approximately 40 throat 
diameters) the surface pressures were low and uniformly distributed with 
the sonic nozzle. As the conic nozzle area ratio was increased, at the 
same distance, the surface pressure also increased. However, with an 
area ratio of 25, the maximum surface pressure was only 0.3 percent of the 
chamber pressure. The shape of the pressure curves with decreasing valves 
of decent height was generally symmetrical, with the peak value on the 
jet centerline. A maximum surface pressure of 5 percent chamber pres- 
sure was recorded for the 25/1 conical nozzle at a descent height of 10 
throat diameters. This nozzle had a throat diameter of 0.50 inch, an 
exit diameter up to 2.5 inch (expan. ratio =s 25/1) and a maximum length 
from throat to exit of 3.75 inches. 

b. Surface pressure distributions were also dependent on 
nozzle contour. Short bell-shaped nozzles gave annular surface pressure 
distributions as a result of shock coalescence originating within the 
nozzle . 


c. The high pressure area on the simulated lunar surface was 
small in comparison with the diameter of the billowing jet. In general, 
the high pressures were contained within a circle whose diameter was about 
16 throat diameters, regardless of area ratio or nozzle contour. 
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d. The surface pressures increased rapidly as the vehicle ap- 
proached the simulated lunar surface. The maximum pressure at a descent 
height of 40 throat diameters was only 0.4 percent of the chamber pressure, 
but increased to 6 percent at 13 throat diameters. 

Data from the above series of tests and hopefully, some future tests 
of a similar nature relating surface pressure to nozzle height above the 
surface, chamber pressure, and nozzle parameters will be utilized in 
Phase 1 of follow-on studies to establish equations relating these vari- 
ables to the ’’effective" momentum energy available in the gas stream 
upon striking the lunar surface. Of course, this "effective" value will 
be less than the theoretical valve because of jet expansion effects and 
loss of kinetic energy in passage of the gas stream through the standoff 
shock wave. 

In addition to the tests described above, Lewis Research Center 
conducted qualitative investigations in a 4 x 6 foot vacuum facility to 
visually explore the erosion caused by the impingement of an under-expanded 
jet on a thick layer of 0.010-inch-diameter, white-sand particles. "Cold 
air at a pressure of 150 pounds per square inch absolute was supplied 
to a small nozzle having a throat diameter of 1/32 inch and a correspond- 
ing weight flow of 0.00224 pound per second. A conical nozzle with an 
area ratio of 12 and contoured nozzles with area ratios of 150 and 500 
were tested at descent heights ranging from about 200 to 600 throat 
diameters. The experimental runs lasted for 2.7 seconds, during which 
time the vacuum tank pressure increased from an initial value of 8 x 10”^ 
millimeter to 8 x 10" 1 millimeter of mercury. Nozzle total-pressure 
ratio thus dropped from 100,000 to 10,000 during this cycle. Erosion 
patterns were recorded photographically through a glass window on one 
end of the tank". 

"During the initial start of the flow from the nozzle a curious 
erosion pattern was obtained. The jet eroded an annular depression in 
the sand, while the central area remained intact. Surface erosion patterns 
of this type have also been observed from the downwash of VTOL aircraft. 

The erosion progressed toward the center with time; a central peak formed 
and projected briefly above the surface and then collapsed. This phenomenon, 
which generally occurred during the first 0.5 second, was more pronounced 
at the lower descent heights where the interior peak was higher and lasted 
longer. The jet continued to "dig" a large circular concave hole after 
the elimination of the central peak. 

The final hole depth and diameter were recorded photographically, 
and the following general results were obtained: 

1. Variations in nozzle area ratio from 12 to 500 had little measur- 
able effect on either the depth or diameter of the hole. 
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2. As descent height decreased from 575 to 192 throat diameters, 
the hole diameter increased from 200 to 250 throat diameters while the 
depth remained essentially constant at about 60 throat diameters. 

3. The sand particles were thrown upward and outward in a sheet 
forming essentially a conical surface with a semivertex angle that aver- 
aged about 50 degrees. However, since the particles apparently followed 

a ballistic trajectory, at no time did any of the sand strike the vehicle. 
Maximum particle heights were estimated to be about 480 throat diameters. n 

Result 1 above can be justified by taking a look at the test nozzle 
characteristics : 



W = 0.00224 lb/sec 

As given in Ref. 14, Equation 4.17, 



where A* = throat area 

P = 150 psia 
o r 

W = 0.00224 
K = 1.4 for air 
R =53.3 ft lb f/lb °R 

t 

o 


745° R 
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at a Mach number of unity, 

T = 0.8333 (745) = 620° R 
P = 0.528 (150) = 79.3 psia 

_ 150 (144) _ 3 

p o " 53.3 (745) “ 0-545 lb/f 

Pl = 0.634 ( p ) = 0.345 lb/ft 3 

C = speed of sound at T = 7 KRT = 1340 ft/sec 


c It 

~ = J = 7 0.832 = 0.912 
o V o 


= 0.912 (C 0 ) = 1220 ft/sec velocity at Throat 


using tables for compressible isentropic flow, the following was obtained 
For A 2 /A = 12 (Case 1) 

P 2 = 0.82 psia 
M = Mach No. = 4.133 
T 2 = 0.2258 (T q ) = 168° R 
v 2 

M* = 2.154 = 35 = Cj_ 


V 2 = 2.154 (1220 ft/sec) = 2620 ft/sec 
For A /A = 500 (Case 2) 

P 2 = 0.004155 psia 
M = 9.828 
T 2 = 36.8° R 
M* = 2.388 

V 2 = 2.388 (1220) = 2920 ft/sec 
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Therefore, although the exit Mach numbers for the two cases vary 
considerably, the exit velocities change a very small amount. Thus, 
there is very little change in the kinetic energy of the exhaust gases 
and as such, nozzle area ratio variations would not be expected to cause 
much change in the "cratering” effect. 

Result 2 indicates an increase in hole dia. from 6.25 to 7.8 inches 
when the nozzle was lowered from 18 to 6 inches. Diameter increase is 
attributable to increase surface pressure with a corresponding increase 
in kinetic energy transmitted/unit surface area; all as a result of re- 
duced jet expansion effects at reduced heights above the surface. Very 
little change in depth may have been a result of compaction of the sub- 
surface particles after initial impingement and also because of the in- 
crease in the test chamber pressure from 0.08 to 0.8 mm Hg during the 
2.7-second test run. Such an ambient pressure increase could have materi- 
ally affected the degree of energy transfer to the surface particles. 
Actually, this test data has practically no quantitative value because 
of lack of surface pressure distribution information. Also, test nozzles 
of the size used are very sensitive to flow variables and can very easily 
give erroneous or misleading test data in addition to being very difficult 
to obtain scaling factors for. It is apparent that much more test data 
taking into account all the dependent variables is required to enable 
future determination of particle motion parameters and cratering effects. 
This must be done by holding all variables constant but one until all 
combinations have been analysed. 

Result 3 indicates that there will probably be no hazards to the 
vehicle as far as backflow of particles is concerned. A maximum particle 
height (0.01 inch dia. particle) of 480 throat diameters, assuming a 
departure angle of 50 degrees from the vertical (following the sides of 
the crater) is equivalent to an initial velocity of 14 ft/sec. Again, 
there is insufficient data to attempt a correlation of these results to 
other nozzle configurations. More testing is in order. 
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FIGURE 6. MOTION OF A PARTICLE IN A CENTRAL FORCE FIELD OF INVERSE SQUARE INTENSITY 
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FIGURE 7 INERTIAL TRAJECTORIES IN A CENTRAL GRAVITY FIELD 
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